Strong pinning and vortex energy distributions in single crystalline Ba(Fei_a^Coa;)2As2 
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The interrelation between heterogeneity and flux pinning is studied in Ba(Fei-xCocc)2As2 single 
crystals with widely varying Co-content x. Magnetic Bitter decoration of the superconducting vortex 
ensemble in crystals with x = 0.075 and x = 0.1 reveal highly disordered vortex structures. The 
width of the Meissner belt observed at the edges of the crystals, and above the surface steps formed 
by cleaving, as well as the width of the intervortex distance distribution, indicate that the observed 
vortex ensemble is established at a temperature just below the critical temperature Tc. The vortex 
interaction energy and pinning force distributions extracted from the images strongly suggest that 
the vortex lattice disorder is attributable to strong pinning due to spatial fluctuations of Tc and 
of the superfluid density. Correlating the results with the critical current density yields a typical 
length scale of the relevant disorder of 40 - 60 nm. 



PACS numbers: 

I. INTRODUCTION 

Recent vortex imaging studies of iron-based supercon- 
ductors have unveiled highly disordered vortex structures 
that challenge the use of traditional analysis procedures 
based on the characterization of positional and orienta- 
tional lattice correlations .^^^ For example, the combina- 
tion of small-angle neutron scattering experiments with 
Bitter decoration^ ^ and Magnetic Force Microscop}^^ re- 
vealed a "vitreous" phase in Ba(Fei_a;Coa;)2As2 single- 
crystals. The latter work shows that the vortex struc- 
ture of the overdoped material {x = 0.19) presents, at 
best, short-range hexagonal order in the field range of 
10~^ to 9T. Disordered vortex structures were also ob- 
served by means of Bitter decoration in single-crystals 
of other iron-based pnictide superconductors including 
Bai_a;Ka^Fe2As2, Sri_2^Ka^Fe2As2, and SmAsGi-a^F^^.^ 
Regardless of material, doping, and synthesis method, 
the disordered vortex structures are attributed to a 
strong pinning the nature of which was not discussed. 

The only reported ordering effect on the orientation of 
the vortex ensemble is that induced by twin-boundaries 
in Ba(Feo.949Coo.o5i)2As2.^ This scanning Superconduct- 
ing Quantum Interference Device (scanning SQUID) mi- 
croscopy study shows that vortices avoid twin bound- 
aries acting as a barrier for vortex motion.^ These re- 
sults echo earlier work on the doping-dependence of the 
critical current density jc, that suggests that structural 
domain walls may act as effective pinning centers in the 
underdoped material.^ 

Concerning the nature of the strong pinning ubiqui- 



tous to iron-based superconductors, inhomogeneities in 
the dopant ions distribution was suggested to be at the 
origin of a dense vortex pinning nanostructure in the 
case of Ba(Fei_a^Coaj)2As2 {x = 0.1).^ The same study 
shows that since thermal fluctuations are weak, the finite 
width of the superconducting transition can only be as- 
cribed to an inhomogeneous Tc distribution due to local 
compositional variations.^ Furthermore, scanning tun- 
nelling spectroscopy studies in different iron-based pnic- 
tides reveal nanoscale variations of the local supercon- 
ducting gap.l^^HI^In Ba(Fei_ajCoa^)2As2, the length scale 
on which the deviations from the average gap value occur 
is comparable to the average distance between dopant 
ions.^^ Nevertheless, no correlation between the vortex 
positions and the superconducting-gap inhomogeneities 
or other defects has, as yet, been found. ^ 

Hence, all techniques agree on the absence of an or- 
dered vortex structure in iron-based superconductors. 
However, there is no clear consensus on the origin of the 
disorder in the vortex ensemble and the pinning causing 
it. The aim of this paper is the characterization of this 
strong pinning by means of a quantitative analysis of the 
spatial distributions of pinning energy and pinning force. 
We found that the key to understand the disordered vor- 
tex configurations is that these are frozen at T ~ Tc, in 
crystals with important spatial variations of the super- 
conducting parameters. The correlation of the extracted 
pinning forces and energies with measurements of jc in- 
dicates that spatial variations of the superfluid density 
and of Tc , on the scale of several to several dozen nm, 
are the most relevant for pinning. 
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FIG. 1: (Color online) (a) Photograph of 

Ba(Feo.925Coo.o75)2As2 crystal # 2. (b) Scanning elec- 
tron micrograph of the decorated sample #2.1 cut from 
the larger crystal #2. (c) Magneto-optical images of 
Ba(Feo.925Coo. 075)2 As2 crystal #2 at T = 15 K and indi- 
cated values of the applied magnetic field, (d) Differential 
Magneto-Optical (DM0) images in the vicinity of Tc for 
fioAHa = 0.1 mT. (e) DM0 image of Ba(Feo.9Coo.i)2As2 
crystal #1 at full screening (T = 17.5 K) and (f) at 
mid-transition (T = 19.25 K). The arrows indicate regions 
of paramagnetic transmittivity at the superconducting 
transition, (g) Local transmittivity Th measured on the 
three regions of Ba(Feo.925Coo.o75)2As2 crystal # 2. (h) Th 
measured on the three regions of Ba(Feo.9Coo.i)2As2 crystal 
#1 indicated in (e). Scale bars correspond to a length of 100 
/xm unless indicated otherwise. 



II. EXPERIMENTAL DETAILS 

Single-crystals of Ba(Fei_a^Coa^)2As2 were grown using 
the self- flux method. Starting reagents of high-purity 
Ba, FeAs and Co As were mixed in the molar ratio 1:(4- 
x) :x, loaded in alumina crucibles and then sealed in evac- 



uated quartz tubes. For each doping level, chemical anal- 
ysis by an electron probe was performed on several crys- 
tals yielding the Co content within 0.5% absolute accu- 
racy. For this work we studied six doping levels. 

The penetration of magnetic flux into selected crys- 
tals of thickness 30 /im was visualized by the magneto- 
optical imaging (MOT) method. A ferrimagnetic gar- 
net indicator film with in-plane anisotropy is placed on 
top of the sample and a polarized light microscope is 
used to observe it. The Faraday rotation of the indi- 
cator allows the detection of regions with non-zero per- 
pendicular component of the magnetic flux density B±^ 
revealed as bright when observed through an analyzing 
polarizer. Dark regions correspond to B± ^ 0. In order 
to characterize the inhomogeneity of the crystals in the 
vicinity of the critical temperature we use the differential 
magneto-optical method (DM0). Images acquired at 
applied fields Ha + Ai^^ and Ha are subtracted, and the 
differential images averaged by repeating the procedure 
50 times. In the present experiments jj^oAHa = 0.1 mT 
(with /io = 47r X 10"'^ Hm"^). 

The field dependence of the global critical current den- 
sity of selected crystals was obtained from magnetization- 
loop measurements conducted using a Quantum Design 
SQUID magnetometer. The critical current densities 
were extracted using the Bean-critical state model. As 
discussed below, the assumption of this model is justi- 
fied by the way flux penetrates into the crystals. Within 
the Bean model, jc = 3A4/Va^ where Ai is the magnetic 
moment, V is the sample volume, and 2a the sample 
width.ISl 

For the Bitter decoration experiments,^^ rectangles of 
dimension 200 /im x 300 /im were cut from larger crystals 
using a 20 /im wire saw and 1 jam SiC grit. Bitter deco- 
rations were only performed on crystals with x = 0.055, 
X = 0.075, and x = 0.1. The sample surfaces were freshly 
cleaved before the experiments (Fig. [l^, b). The exper- 
iments were carried out at liquid Helium temperature 
(4.2 K) and He-exchange gas at pressures of the order 
of 200mTorr. The images shown here are the result of 
field-cooling experiments at a field fioHa = 1 mT applied 
parallel to the c-axis of the crystals. The decorated vor- 
tex arrangement was observed by scanning electron mi- 
croscopy at room-temperature. 



III. RESULTS 

A. Magneto-Optical Imaging and jc measurements 

Figure [l](c) shows examples of magneto-optical im- 
ages, here obtained at T = 15 K on single-crystal #2 of 
the composition with x = 0.075. The images reveal a 
globally homogeneous penetration of the magnetic flux 
into the sample obeying the Bean critical state. ^^'^ We 
obtain the local value of the critical current density from 
jc ~ 6 dB± I dx (the factor 6 is estimated from Ref . |14] 
for a crystal aspect-ratio of 0.1). The DM0 images in 
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FIG. 2: (Color online) (a) Transition temperature, Tc, versus 
Co doping-level. The error bars denote the local spread of Tc 
values within a given crystal. For each doping level, #1, #2 
and #3 denote different crystals. For x = 0.075 and 0.1, 
the numbering denotes that of the decorated crystals, (b) 
Co doping-level dependence of the critical current density jc 
measured by MOI at B = 30 mT and a reduced temperature 
of T/Tc = 0.47. 



reported in the literature for the same material P 



B. Vortex imaging 

The Bitter decoration technique^^ was used to observe 
vortex structures on three of the crystals used to compile 
Fig. [2j more precisely, on crystal #1 of the composition 
with X = 0.1, crystal ^2 with x = 0.075, and on crystal 
7^2 with X — 0.055. The decoration of crystal with 
X = 0.055 was unsuccessful, presumably due to the large 
value of the penetration depth at low doping. The deco- 
rated patterns reveal highly-disordered vortex structures 
as in Refs. |l|3f [6). Figures |3](a) and (b) reveal regions 
of high and low vortex density, as well as the formation 
of vortex-free zones near the crystals edges and surface 
steps, due to the circulating Meissner current. These im- 
ages are representative of those obtained on other regions 
of the crystal surfaces after different cleavage runs, and 
on other crystals. From the images, we extract the av- 
erage value of the magnetic induction as Bint = ^v^o, 
where riy is the vortex density and = h/2e is the flux 
quantum. For all images we obtain an average induc- 
tion Bint ^ 0.8 mT, 20% smaller than the applied field 
Ha = lmT. 

Figures [4](a) and (b) present the Delaunay triangula- 
tions of the images in Fig. |3]for x = 0.075 and x = 0.1, 
respectively. Here, the blue dots represent vortices with 
sixfold coordination while the red dots represent vortices 
which have a different coordination number. The insets 
to Figs. [4](a) and (b) show the Fourier transforms of the 
vortex positions which once again demonstrate the ab- 
sence of any order in the vortex structure. 



Fig. [T](d) reveal the same Bean-like flux penetration with 
an inhomogeneous jc arising from the spatial variation of 

This inhomogeneity can be quantified using a plot 
of the local transmittivity, defined as the ratio Th = 
[I{r, T) - /(r, T « T,)]/[I{v, T » T,) - I{r, T « T,)] of 
the relative local luminous intensities /(r, T) in the DM0 
images. The temperature-dependence of Th measured 
on different regions of crystals #2 and #1 is depicted in 
Figs . [l](e),(f). The local variation of Tc-values within a 
given crystal is of the order of 0.5 - 1 K. In addition, re- 
gions of lower Tc give rise to a paramagnetic signal at the 
transition due to flux concentration by the surrounding 
superconducting parts of the crystal. 

Figure |2](a) summarizes the width of the Tc distri- 
bution for a large number of Ba(Fei_a^Coa;)2As2 single- 
crystals of different doping levels. Figure |2](b) shows the 
Co doping-level dependence of jc for the same series of 
single-crystals at a reduced temperature of T/Tc = 0.47. 
A rather large sample-to-sample variation of the low-field 
{B_\_ = 30 mT) jc is observed. Certainly, no clear doping- 
dependent trend appears, as proposed in Ref. |8j. The 
obtained critical-current values are comparable to those 



C. Vortex configurations near surface steps 

The correct determination of the distribution of vortex 
pinning energies in the crystal and its interpretation re- 
quires knowledge of the temperature at which the vortex 
ensemble was frozen in the observed configuration. To 
determine this, we analyze the vortex distribution near 
the ubiquitous steps seen on the surfaces of the crys- 
tals. Such steps result from the repeated crystal cleav- 
age performed during the Bitter decoration experiments. 
In zero-field cooled experiments, steps act as obstacles 
for vortex entry into the sample; they were described in 
Ref. 121] as "vortex diodes". 

However, the present decoration experiments are in 
field-cooled conditions and hence vortices nucleate in the 
sample at the same temperature that the mixed state is 
stable. As one cools down, the Meissner screening cur- 
rent running along the crystal edges, but also along the 
surface steps, increases as the penetration depth Xab for 
currents running in the a6-plane decreases. Thus, while 
cooling, vortices on the high-side of the step are progres- 
sively repelled by the increasing Meissner current density 
jM ^ Ha/ Xab- At the same time, the proximity of the 
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FIG. 3: Bitter decoration images of Ba(Fei-xCox)2As2 single- 
crystals (a) crystal #2.1 with x = 0.075, and (b) crystal #1 
with X = 0.1. The graph above panel (a) shows a profilometer 
measurement upon crossing the step that traverses the image 
from top to bottom; full vertical scale is 2.5 /im. The width u 
of the Meissner belt behind the step is also indicated in (a). 



step surface results in an attractive force that can be de- 
scribed by an image vortex segment. Finally, the vortex 
lattice elasticity tends to restore a homogeneous flux dis- 
tribution near the step. The situation is therefore similar 
to vortex entry or exit over a surface barrier. 

At the low fields of interest, the single vortex part of 
the tilt modulus dominates vortex elasticity, so that the 
force balance is 
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IS the vortex line tension, = 
IS the vortex line energy, sx = Xab/^c is 
the penetration depth anisotropy, v = u/Xab^ represents 
the width of the vortex- free zone behind the step, nor- 
malized to Xab' The step running through Fig.[3](a) has 




FIG. 4: (Color online) Delaunay triangulation of the vortex 
structure of Ba(Fei-xCocc)2As2 single-crystals (a)# 2.1 with 
X = 0.075 and (b) #2 with x — 0.1. The blue dots represent 
vortices with 6 nearest neighbors while red ones represent vor- 
tices with a different coordination number. The insets show 
the respective Fourier transforms of the vortex positions. 



a height h = 1.5 /im while the vortex-free region behind 
it has a width u = 1.8 /im. Estimating the penetration 
depth anisotropy £x ^ 0.16 from Refs. i23l and 24, and 
with all other parameters known, Eq. [l] can be solved 
graphically to yield v ~ 1.5, that is, A^^ ^ 0.6ii ~ 1.2 jim. 
Combining A^b-data from Refs. [6| and |23l we conclude 
that the observed vortex pattern is frozen at Tf ^ 0.9Tc. 




FIG. 5: Representation of vortex lines near a surface step 
under field-cooled conditions 
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FIG. 6: (Color online) Left-hand panels: Normalized color- 
coded maps of the vortex interaction energy calculated from 
the images of Fig. |3]in Ba(Fei-ccCoa;)2As2 single-crystals with 
(a) X = 0.075, and (b) x = 0.1. Right-hand panels: His- 
tograms of the normalized interaction-energy histograms for 
(c) X = 0.075, and (d) x = 0.1. 
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FIG. 7: (Color online) Color-coded maps of the modulus of 
the individual vortex pinning force per unit length, calcu- 
lated from the images in Fig. [sjfor Ba(Fei-xCox)2As2 single- 
crystals with (a) X = 0.075 and (b) x = 0.1. (c) and (d) 
represent the pinning force distribution for x = 0.075 and 
X = 0.1 respectively. 



D. Pinning energies 

The inter- vortex interaction energy is calculated from 
the vortex positions obtained from the decoration images. 
We calculate the interaction energy 



(2) 



per unit length along the vortices' direction. Ko{x) is the 
lowest-order modified Bessel function, and the vortex line 
energy sq oc is proportional to the superfiuid density. 
We take into account all vortices j situated at a distance 
smaller than 10 Xab from vortex i. This cutoff radius was 
chosen after verifying that the interaction energy does 
not change significantly if greater values of j are consid- 
ered. For the determination of the energy distribution 
histograms, we only take into account vortices situated 
away from the edges of images, at a distance larger than 
4Aa6.^^ We used in this calculation the penetration-depth 
value at the temperature at which the vortex structure 
was frozen, KhiT /T^ = 0.9) (see Section |IIIC ). 

A similar procedure yields maps of the pinning force 
acting on an individual i vortex per unit length 
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with Ki{x) the first-order modified Bessel function. 
Since the system is stationary, Newton's third law re- 
quires the repulsive force exerted by neighbor vortices be 



balanced by the pinning force. A map of the modulus 
|f^| thus represents a map of the minimum pinning force 
acting on each vortex. In the case of a perfect lattice 
resulting from negligible effect of pinning, the sum (|3| 
vanishes. 

We present our results by color-coded maps spanning 
the whole decoration images of Figs.[3](a) and (b), and 
by histograms of the interaction energy distribution. The 
interaction energy maps with the energy-scale normalized 
by So are shown in Figs. [6](a) and (b). A granular struc- 
ture of denser regions with larger interaction energy, and 
dilute regions with smaller Sint is clearly visible. This 
granularity is translated into broad vortex interaction- 
energy histograms as shown in Figs. [6](c) and (d). The 
histograms are reasonably well fitted by a Gaussian dis- 
tribution. The standard deviations of these histograms 
are of the order of 23 %, in contrast with 50 % for the 
rather regular vortex structured of the same density 
imaged in the high- Tc material Bi2Sr2CaCu208+5. How- 
ever, as a result of the high reduced temperature Tf /Tc 
at which the vortex ensemble is frozen, the mean inter- 
action energy (normalized by Sq) is ten times larger in 
Ba(Fei_a,Co^)2As2 than in Bi2Sr2CaCu2 08+5!^ 

The reduced temperature Tf/Tc at which the vor- 
tex ensemble is frozen not only affects the deduced in- 
teraction energies, but also has a profound effect on 
the (orientational) order observed in the decorated vor- 
tex ensemble.^- Pardo et al. reported^^ that in op- 
timally doped Tl2Ba2Cu06-5 superconductors with a 
broad magnetically reversible regime in the temperature- 
field phase diagram, and concomitantly low Tj/Tc, Bitter 
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decoration yields a regular triangular lattice, while dec- 
orated vortex ensembles in the overdoped material with 
a narrow reversible temperature range (and high Tj/Tc) 
are amorphous. At the origin of this effect is the high mo- 
bility of vortices just above Tf in materials with a wide 
reversible regime, such as Bi2Sr2CaCu208+5 or optimally 
doped Tl2Ba2Cu06-5. On the other hand, the low mo- 
bility of the vortices just above Tf due to strong pinning 
in the vortex liquid phase in materials [such as, appar- 
ently, Ba(Fei_a^Coa;)2As2] that have a narrow reversible 
regime yields an amorphous vortex ensemble. 

Fig. [7| shows maps of the modulus of the pinning 
force of individual vortices per unit length of the vor- 
tex lines. The pinning energy shows some correlation 
with the interaction energy at the local scale: regions of 
large (small) Eint generally correspond to regions of large 
(small) |fi|. There is noticeable inhomogeneity on scales 
smaller than the apparent grain size. The juxtaposition 
of a region with homogeneous large superfluid density 
{i.e. So) with a region of homogeneous small Eq would 
give rise to a larger pinning force at the interface only. 
In the images, fluctuations of the pinning force within 
grains of similar Sint are observable. Therefore, inho- 
mogeneity of the superconducting parameters exists not 
only on the /im scale of the images, but also on smaller 
length scales. 

It is interesting to note that the rendered pinning forces 
are simply related to a met ast able current density j^, run- 
ning through each vortex, as = (<l>o/|B|)B x j^. The 
average pinning force per unit length of 5 x 10~^N/m, 
with local maxima of up to 6 x 10~^N/m, imply local 
currents of the order of 2.5 x 10^ Am~^. Maximum cur- 
rents are of the order 3 x 10^^ Am~^, comparable to the 
low-temperature value of the critical current density. 



IV. DISCUSSION 

Since the vortex locations result from the balance be- 
tween inter- vortex repulsion and the interaction £p of in- 
dividual vortices with the pinning impurities, one has, 
at Tj, Sint = Ep- The position of the maximum and 
the width of the interaction-energy distributions [see 
Figs. |6](b) and (d)] are therefore determined by, respec- 
tively, the mean and the standard deviation of the pin- 
ning energies, at Tj, of the individual vortices in a given 
image. In particular, the displacement of the maximum 
of the distribution with respect to the position of the S- 
peak energy-distribution of a perfect vortex lattice of the 
same density is a measure of the mean pinning energy. 
As far as the vortex densities of Fig. [6] are concerned, 
the average Bint = 0.8 mT yields a (5-peak- maximum at 
Sint = 2.560 • By comparison, the maxima of the distri- 
butions for both investigated crystals in Figs. [6](a) and 
(b) occur at Sint ^ 3.2eo. The average pinning energy 
per unit length is therefore Sp ^ O.Tsq, while the variance 
in pinning energy is given by the width of the distribu- 
tion, {{S^) - {Spfy^'^ - 0.5£:o. Note that 3.260 corre- 



sponds to the interaction energy of a triangular vortex 
lattice with ^on^ = 1 mT, i.e. the external field applied 
during the experiments. This means that the average in- 
teraction energy is determined by vortex-rich areas, with 
<l>on^ ^ 1 mT. However, the vortex density also presents 
vortex-poor areas so that the average Bint = 0.8 mT. 

The large absolute values of the inferred pinning en- 
ergies can be understood if one combines the notion 
that the crystals show local variations both of the crit- 
ical temperature Tc = T'c(r) and of the line energy 
^0 = ^o(r,T), and that Tf /Tc ^ 1. As the crystal is 
cooled below Tc, vortices will avoid regions of higher 
Tc and ^q, and accumulate in regions with lower values 
of these parameters. They will remain trapped in such 
regions as the temperature is lowered below Tf. The 
large absolute values and variances of the pinning ener- 
gies revealed by the decoration experiment are caused by 
the local variations of Tc(r), which manifest themselves 
through the temperature dependence of the line energy, 
6:o(r,T) = £o(r,0)[l - T/Tc(r)]. More specifically, the 
width of the inferred pinning-energy distribution (Fig. [6| 
should correspond to the width Aso{Tf) of the line en- 
ergy distribution, 

0.5eo(0)(l - Tf/T,) ~ A£o(I». (4) 

Near to the critical temperature, Aeo{T) = 
eo{0)TATc/T^ is determined mainly by the width 
ATc of the distribution of local Tc(r). Solving Eq. Q 
then yields Tf = Tc/[1 + ATc/0.5Tc]. Taking Tc = 24 
K, and estimating ATc ^ 0.8 K from the DM0 data of 
Fig. [Tj one obtains a freezing temperature Tf = 0.94Tc 
for X = 0.075; the same exercise with Tc = 19 K and 
^Tc = 0.5 K yields and Tf = 0.95Tc for the crystal with 
X = 0.1. Thus, the analysis of the inhomogeneous and 
disordered vortex distribution, as well as the vortex 
distribution near steps and edges, is fully consistent 
with the observed patterns having been frozen between 
T = 0.9 and 0.95Tc. We can draw the same conclusion 
from the local variations of the vortex density. For 
example, for the crystal #2 with x = 0.075, the largest 
local vortex gradient correspond to 0.15mT//im or 
1 X 10^ Am~^. This value is consistent with the critical 
current density of Ba(Feo.925Coo.o75)2As2 crystal #2 at 
23 K. 

At low temperatures, the spatial variations of the mag- 
nitude of the line energy £o(r,0) are dominant. These 
correspond to the variations in space of the superfluid 
density,^ and are responsible for the non-zero low-T pin- 
ning force associated from spatial inhomogeneity. A spa- 
tially homogeneous superfluid density would imply a van- 
ishing (or logarithmically weak) pinning energy at low T, 
at odds with the existence of a large critical current den- 
sity (see, e.g. , Fig. [8|. As in all charge-doped single- 
crystalline iron-based superconductors, the critical cur- 
rent of Ba(Fei_a^Coaj)2As2 is composed of a contribution 
from strong, extrinsic pins, and from a contribution from 
pinning by atomic sized-point pins. The latter dominates 
at high fields [here, above IT at 5K, and above 0.2 T at 
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FIG. 8: (Color online) Critical-current densities in 
our Ba(Feo.925Coo.o75)2As2 crystals. (a) Temperature- 
dependence of the low- field jc of crystals #1 (x = 0.1) 
and #2.1 (x — 0.075), as obtained from MOI. Error bars 
represent the dispersion of jc within a given crystal, (b) 
Field-dependence of jc for crystal #2.1, obtained from mag- 
netic hysteresis measurements using a SQUID magnetometer. 
Straight lines indicate fits with Eq. ([6|, see section IV 



17.5 K, see Fig. [8[b)],^^ while the former contribution 
manifests itself as a low-field plateaJl^ 
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The crossover field 5* is that above which the num- 
ber of effective pins per vortex is limited by the in- 
tervortex repulsion, fp is the maximum pinning force 
exerted by a single strong pin, rii is the pin density, 
and /7p/[J] is the pinning energy gained by a vortex 
line traversing such a pin. The measurement of the 
low-field critical current density jc(0) and the slope 
djc{B)/dB~^^'^ allows one to eliminate rii and to obtain 

fp = TT^l^^ex {jc{^)/[djc{B)/dB-^/^]} from experimen- 
tal data without further assumptions. We find that, at 
5 K, ^ 3 X 10-13 N for both Ba(Feo.925Coo.o75)2As2 
crystal #2.1 and Ba(Feo.9Coo.i)2As2 crystal #1. 

The identification of the strong pins with regions of 
lower eo{T) means that fp should be interpreted in 
terms of the local maxima of the position-dependent 
force /(r) = J^^ \/£o{'^^)dz experienced by vortices as they 
move through the sample. Here Sz is the maximum ex- 
tent of a region of low eo{T) along the field direction. We 
approximate 



(7) 



where Sg is the length scale characterizing the disorder 
in the direction perpendicular to the field, and A^o is 



the standard deviation of the eo{r) distribution in the 
crystal. The pinning energy Up ^ fpSg. A comparison 
of Eq. ^ with the value of fp obtained from jc yields 
A^o ~ 3 X IQ-i^ Jm-i for a unit aspect ratio 5z/5g. 

In a next step, we evaluate the ratio of fp/\^i \ to obtain 
the average distance between effective pins, £ = 60 nm. 
Using Eq. (17) of Ref.l30l which has C = [sxliiniUpfl^ , 
one finds (nibz)~^l^ ^ 60 nm. With all parameters 
known, the low-field value of the critical current density 
is reproduced as 



Jc 



rl/2 



AgQ 



Sg\l eo 



(8) 



in fair agreement with the data of Fig. [8](a). The inves- 
tigated features of vortex pinning in Ba(Fei_a^Coaj)2As2, 
including the disordered vortex patterns and the critical 
current density, are therefore consistently described by 
the presence of spatial variations of the superfluid density 
on the scale of several dozen nanometers, in agreement 
with the conjecture of Ref. [9j. 

Note that the observed spatial structures at the macro- 
scopic (Fig. [T]) and mesoscopic (Fig. [3| levels are not 
those responsible for the critical current. The random 
vortex positions observed in the decoration experiments 
are determined by the underlying nanoscale disorder, an 
observation consistent with the fact that disordered vor- 
tex structures have been observed up to high fields.'^ 

One may speculate about the possible link between 
the existence of nm-sized regions of reduced superfluid 
density, the local variation of the dopant atom density, 
and the effect of the overall doping level. For instance, 
one would expect the fluctuations of the Co density to 
be more important at lower doping levels, yielding larger 
local fluctuations of Eq. However, given the much larger 
values of the penetration depth at low doping, we have 
not been successful in performing Bitter decorations on 
the relevant crystals. Recent STS studies have reported 
substantial variations of the v alue o f the superconducting 
gap on a scale of 10 to 20 nm.^ These local variations 
of the gap magnitude should correspond to the variations 
of the critical temperature and therefore lead to vortex 
pinning. Although it is tempting to relate our results to 
the nanoscale disorder observed in the STS gap-maps, 
it should be remarked that the spatial scale of the 
variations in the gap maps is a factor of 3-6 smaller than 
that found from the analysis of the data presented here. 
This would correspond to a concomitantly larger jc in 
the samples used in Refs. f 7|ll| . 



V. CONCLUSION 

Bitter-decoration imaging of the disordered vortex dis- 
tribution in superconducting Ba(Fei_ajCoaj)2As2 single- 
crystals with X = 0.075 and x = 0.1 reveals a substan- 
tial local variation of pinning energies and pinning forces. 
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The magnitude of these fluctuations is suggested to stem 
from nanoscale spatial variations of Tc and/or the su- 
perfluid density due to an inhomogeneous distribution of 
dopant atoms. The spatial scale of the variations is in- 
ferred from the correlation of the features of the vortex 
distributions with global and local critical current den- 
sity measurements. The macroscopic spatial variations of 
the critical temperature observed using magneto-optical 
imaging give an idea of the magnitude of the Tc variations 
in the crystals, but are unrelated to the measured pinning 
properties. The same can be said for mesoscopic disorder 
structures observed by single- vortex imaging. An impor- 
tant corollary of our work is the fact that the observed 
vortex distributions are frozen, at a length scale of the 
lattice spacing, at a high temperature close to Tc. 
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